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Abstract 


The puiposp ci tt.is (-Sffr is tu suprarize 
results obta)r>rd fe r the acisftivw control of t.>e 
F- ec aircraft jsl' rnt- so- cal ird ^KKAC retnod. 

The discusstcr. <r.cl..es tc.e selection cf the p*r~ 
forear re criteria ft tictr. tre lateral ard ♦.‘■e 
lor.91tviir.al dvrar.es, the design of the Kal.-aui 
filteis for iffeicr.r f;i-j..t crr.ditions , the 
“ider.tif itatitr “ asp ets r » the design using 
hypothesis te*tir.9 liras, ar.d t.-.e fsrforrarce of 
the closed loof sdaptivc syster..^^ 

* . i - V / aA" ' - ' ' ' 

1. Introduction 

We purpose of this paper is to prasen*. 
preliminary results sr. a svudy irhich involves 
the appllcaiiot. cf advanced adaptive control 
techniques tc the design of a stability augmen- 
tation systcp ir teth tne longitudinal and 
lataral dyr.arics of the F-nc aircraft. NAFA nas 
beer using the F-9C aircraft as a test vehicle 
for evaluating different digital-f ly-by-wira 
(Df'BH) contril tec.tr.iques , using the IBM AP-101 
as the airborne cerputer we rerark that tlie 
eventual ir^lerer.tat ion o: the control aloorithpia 
on tlie specific ai.rorne computer has had a major 
impact upon the pnilo'ophy adopted for the design 
of the control systvr in view of the obvious sto- 
rage and leal-tipc crrputatior.al constraints. In 
addition, the dcsi,r. was crucially dependent upon 
the Bcnsors that could be utiliceo in the sense 
that sensors that utilized external aarodynamic 
awaauresients . e.g., airspeed, altitude, angle of 
attack and sideslip vanes should not be employed 
in the candidate design. Thus , the design guide- 
lines required that the sensors esrociated with 
the adaptive control system shoul'i l<e liauted to 
accelerometers, rate gyros, and perheos attitude 
sensors (although the latter wen deemed undesir- 
able in view of their errors whei the aircraft 
underwent severe pilot induced awncuvers) . 
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■cm the viewpoint of aedelling. It is ob- 
vious the the dynamic state aquations of an air- 
craft involve nonlinear differential aquatione 
(see Etxin (11). Howi-ver, the Inforawtion given 
by NASA Langley Peseerv-h Canter (LRC) to the KIT/ 
ESL team consisted in the specification of the 
uncoupled , linear time-invariant open-loop 
longitudinal and lateral dyi amics of the F-8C 
aircraft associated with eqiiilihriim. flight . 

Table I gives a list of th.c flight conditions 
that were available for the design, ihus. the 
general structure of the <-ouationa were of the 
form x(t)»Ax(t)*Bu(t) . The numerical values 
of the elcsmnts of the A and B matrices can be 
found ' ' a report by Gera |2] , baaed upon wind 
tunnel tests, and a report hy Hooley and Fvan^ 

(3), baaed upon linearization of the nonline.ir 
dynamics employed by NASA/UC for their nonlinear 
simulation of tha F-BC aircraft. He remark at 
this point chat the nimierlcal values for the A 
and ^ mat rices given in [2] and (31 are not iden- 
tical reflecting the feet that different sources 
ware used to obtain them. ITie design reportad in 
this paper is based upon Gera's report [2]. 

The feet that the 16 flight conditions span 
an estreswly wide envelope for operating the 
aircraft, with drastic changes in the open-loop 
dynamics, makes the fixed-gain design of the con- 
trol system unrealistic. Furtherscrc, handling 
qualities requiresients , such as the C* criterion, 
indicate that pilots desire different closed-loop 
dynamics at different flight conditions. Thus, 
some sort of 'adaptive* gain-scheduling control 
aystem was required, itowever , straight-forward 
gain schedull.ig based upon quantities such as 
velocity, altitude, end dynamic pressure was not 
permitted in view of the sensor rostrictlona aien- 
tloned above. Hence, tlie adaptive control system 
had to « designed in a novel way. 

An additional restriction on the design was 
that the sensor noise and wind disturbanc-s had 
to be incorporated. This led to the need for 
eagiloying Kalman filters, with constant coeffi- 
cients because of the coaputcr amsnry limitations. 

The above problem overview aets the ground 
for tlie specific adaptive control technique which 
we selected to investigate in greet detail. He 
cell the adaptive control technique the Hultiple- 
Hodel-Adeptive-Control (MHAC) method, and we shall 
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dtscttsu IX. in BLir* detail in Section % of thin 
por^^' If ** only >jn« of sowrnl tachniquos 
b«:>«d upon 4*v*l<i|i««;nt<i in nndorn (ontrol thoory 
(»»« t)>« «u(wy nrticl* by Atnons wid Vainiyn 
|4|) and it han itn oriqins in -cnbininq 
hypothrtia-'attinq and ttuchastic control tdoas 
(»<*'■ trfrronc** (4) to IB|) It wan naloctod for 
thi» ktudy bacauar of it» {otantlal pioauk# in 
acadanic caarplra |SI*|b]. and tacaoaa ita aanDry 
and ical'tina corrutatiofial raquiiaoanta could be 
readily asaekaed in view of ita non-iterative 
nature . 

As eaplair^d in norc data. 1 in Section S, 
the rWAC Mthod requirea that a full blowr. steady 
atate Linear -Quadrat ic-''-ausaian (L^.) controller 
be ipplcivnted for each flight condition. Thia 
r.eceasitated the develi^pnent of suitable quad- 
ratic perforaiance criteria for both the longitud- 
inal and the lateral dynanicsi these are 
described for the uont mucus tire case |V| in 
Sections 2 and t, respectively. For leipleiaenta- 
tsen, one nee<*:i a discrete tine Ugc controller 
(10). This is described in Section 4, together 
with the discuaaio.' of sensor errors. TTie KMAC 
alqorithii is described in Section S. The siaula- 
tion results usittq the nonlinear F-8C dynaeucs 
arc described in Section f.. Section 7 presents 
the major conclusion of our studies so far. 

He rerark that in this paper we Jhall only 
focus our attention to the regulation aspects of 
the proMen. i.e., return to equilibrluai flight 
from soaw initial conditions and in the pre- 
sence of stochastic wind disturbances. In our 
study we are considering the proper way of incor- 
porating huran pilot inputs tor both the longi- 
tudinal and lataral case. However, we shall not 
present in this paper any of the approachas and 
pralisupary results for the pilot input case. 


2. Longitudinal Dynamics 
2.1 Introduction 
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In this saction we present an overview of 
the LQC philosophy adopted for designing the reg- 
ulator for tha lonaitudinal dynamics. Attantion 
is given in the development of the quadratic per- 
fomance index end the subsequent model simplifi- 
cation using a short period approxisMtion. Ihe 
Slain concept that we wish to stress is that the 
quadratic performance criteria employed changed 
in a natural way with aach flight condition. Tha 
Burprlsing result was that the short period poles 
of the resultant longitudinal closad-loop system 
were characterized for all flight conditions by 
two constant daispinq rstios, one associated with 
ell subsonic flight conditions and one associated 
with all aupersonic flight conditions. 

2.2 The Longitudinal State Description 

Because of s rate constraint saturation on 
the elevator rate, the control variable selected 
was the tiaie rate of change of the coemi ended 
elevator rate (6^(t)). This was integrated to 
generate the actual coamwnded alavator position 




M 


which was introduced to a first ordar 
nnrvo with a tiaa eoiwitant of 1/12 aeconds to 
generate the ectual devletlom of the elewetor 
6 It! from its trlaned vnlua. Tlie elewetor wee 
them releted to the four “necural* lomqltwdinel 
state variables nanely pitch rats, q(tl (rad/aac) , 
velocity error v(t) (ft/sec) , perturbed angle of 
attack from its trianed value, altllrad), and 
pitch attitude deviation fren ita triened value 
B(t)(rad), 12). In addition, a wind disturbance 
state w(t) wee included (see Appendix A). Thus 
the state vector x(t) for the longitudinal dynasm- 
ies was charactarlsad by seven components 

(2.1) s'(t)-(q(t). v(t), o(t), e(t). 6^(t), 
«^(t), w(tl) 

and tha control variable u(t) was tha c oia n dad 
elevator rate 

(2.2) u(t)-4^^(t) 

This lad to a linaar-tisn Invariant charactarize- 
tion for each flight condition of the form 

(2.1) x(t)-^x(t)*Bu(t)*L^C(t) 

where ^(t) was zero mean white noise, geteratlng 
the wind disturbance and accovntina for randoe 
actuator errors. The cianents of ^ and ^ 
changed with each flight condition while 

(2.4) ^(0 0 0 0 0 1 01' 

2.1 Tha Longitudinal Cost Punctiunal 


In ordar to apply tha standard steady state 
LQC procedure |9] a quadratic performance index 
)ias to be selected. The gonersl structure of the 
index was 

(2.5) J-/?'(t)24x(t)4o'(t)R^u(t)dt 

Note that tha weighting swtrlces had to bi.'! 

different from flight condition to flight condi- 
tion reflecting in a natural way t)iat tha pilot 
wants different )iandling qualitiaa as ths speed 
(and dynaadc pressurs) changes. 


In ths Initial design it was decided that 
should relate the msxianas deviations of 
o pitch attitude, 6max 
o pitch rate, qms x 
o norsial acceleration, s^^amx . 

o maximimi c o —an ded elevator rata, 6^^max 
resulting in the following structure of 
the perforsmnee criterion 


(2.6) J| 
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The normal sccelsrstion 9'*« va* »<>t 

used as a stats variable. However, it is linearly 
related to soa« of the longitudinal state vari- 
ables according to ths formula 

(2.7) v(t)+Kia(t)eKi6^(t)j in g*s 


ITAIL is: a ii'i:'-.! 


V* b«in 9 th« •(•••d. Ift* oafi»t«nta 

hl, k(, k| CM b« calrulatvd fra* th« open loor 
A Mtrtcva, Md h*nv* ct>M 4 * with flight comdi- 
Tton. th« atructur* of th« crttorioa 

(2.6) iai'liM Chat t( at c>0 tha ■taai*\*i valuaa 
of accalat.atton, pitch rata, or pitch attitwda 
occur rad, tiian ana rnmld ba wiUiiy/ to aaturata 
tha alavator rata to rrar>v« theoi. For tha pra- 
lioinary daaiyn tha following nuMarlcal valuaa 
«*ara aalactad (with tha ha Ip of T. Elliott and 
J. Cara) 

(2.1) qnaa-lO^/Vr, 0^_"*pAaan 

,6 >0.41% rad/aac 

acoaa 

vhara an la tha (1,3) alaaiant of tha opan loop 
lanpitadlhal ^ aatria. 

Mouahly apaaklna, this critarion aaana tJtat 
ona la ulllinq to aaturata tha alavator rata 
(0.4JS rad/aac for tha F-SC) if a normal accal- 
aration of 6q waa (alt, or a pitch rata aquiva- 
lant to ICq'a, or a pitch ..rror which if 
tranalatad to ar.qla of attack would alao qanerata 
a 6q norma) accalaratior . 

Iha abova nvaa*rical valuaa wara trMalated 
into the approprl <ta C r.atria (non diagonal poa> 
itive aamidef inita) which char.qad freo fliqnt 
condition to fliqht condition, whila R •^•1/ 
(0.415)* for all flight conditiona. Sneo, tha 
raaulting hQC problcn could 1« aolvad uaing 
availaola coaiputpr aufcrourinaa [11]. 

2.4 koducad Longitudinal Design 

Tha design was nodifiad for two reasons. 
First the gain from the velocity state variable 
v(t) was estreamly snail. Second, it was desir- 
able to avoid using tha pitch sensor. The pitch 
0(t) is weakly observable froai the system dynast- 
ies so that even if a Kalman filter was used in 
t)ie a)Mence of pitch neasuresients, large estima- 
tion errors would be obtained which would adveisly 
affect t)ia perfoma. "he control system since 

there it significant .ceuwack from the estisMted 
pitch attitude. At any rate, since a pilot would 
fly the aircraft he would be able to control 
pitch hismelf. 


time tradeoff betwee n the meaimum noeiael eccelera- 
tiom. a„„^, and mesim >» pitch rate, Ihls 

la c eee latent with tha C* criterion (liTT" 

Mhem the ahort-period c loeed- loop polee were 
ewelweted for both deatgna uaing the maerlcel 
values given by (2,d), we foiaid the unespected 
rennlt that the danring ratio was constMt (0.4db) 
for all 11 sidaonic flight conditions , and a lao 
conetant (0.361) for all tha sie wi sonic fllgtt 
cowdtt iona . Tha cloaed-loop naturally fraguancy 
increased with dyneeic pieenure. 

Since no pole-placement techniquea tsere en- 
ployed (i.e., the nethesMtics t^re not told to 
piece the cloaed-loop polee on a conatMt daag>- 
ing ratio llna), we oonetmeted e tradeoff by 
changing (Jecreeeing) the neaisnw pitch rate 

. This would Increase the pitch rate penalty 
In tha coat functional, Md ona would aspect a 
hlghar d aa gi lng ratio. Iha following values of 
were ei^loyed 

(2.10) «2^^-10gA«, BgAi, 6gA», 4gAt 

Onoe nore the constant daag-ing ratio phenonenon 
was observed, i.a., for each value oi <)^w the 
short period closed loop poles for all aissonc 
flight conditions fall on a constant daagiir.q tatic 
line, Md sieilarly for all supersonic flight 
conditions. This was further verified by consid- 
ering M additional 13 different flight condi- 
tions. 

The nuBierical results are presented in Table 
II. The reason for this regularity of t)ia solu- 
tion of t)ia LQ problem is under investigation. 


3. Lateral Dynamics 
3.1 Introduction 

In this section we present the parallel 
philosophy for the developme n t of tha control 
systam for tha lateral dynaadcs. In this case 
the development of a performance criterion was 
not as straight-forward as In tha case of tha 
longitudinal dynasdes. For m axtanslve dlscua- 
aion see the s.N. tliasls by ~reene (131. 


This lad ua to alisdnatlng the velocity error 
v(t) Md pitch 0(t) from t)ie state equatiAis and 
obtaining the "short period" approximation (5 
state variables). Since pitch did not appear tha 
criteria (2.6) was modified to 


(2.9) Jj 


,L0M 




Md the resultMt LOG problem was resolved. 


2.4 


try of Results 


3.2 Tha Lateral Dynasde. Stats Modal 

The control variablea aalectad for lateral 
control were 

(3.1) U|(t)wj (t)»rnsmanded aileron rata 

~ (rad/acc) 

(3.2) 111 (t)»j (t)-cnmii.1sd rudder rate 

(rad/sec) 

so that tha control vector la defined to be 


From the viewpoint of translMt responses to 
t)ic variables of interest (normal acoalaration, 
pitch rate, Mgle of attack) tha transient re- 
sponses to initial conditiona wara alxwst identi- 
cal for both designs. Thus, tha short period 
motion of the aircraft wm dominated by tha rela- 


(3.3) u'(t)-Iui(t) u;(t)J 

The servrxwchanlcs were taken Into account. Tha 
coMmnded aileron and rudder rates' wara integrat- 
ed to generate tha coam ended aileron (6^^(t)) and 
rudder «^g(t)) positions, rsspectivsly. For ths 
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r-tr the cuMMr.drd •ilaron r«t« 

« flrat ordar laa ftvivo, with m ttna con* 
•tritet of l/k> Mw 'AsM, tc 9CfMrrat« th* actual 
atlarwr (oaition {,,(1' (rtd* ) . 

n>« ccwMndad ruddar latr £^f(t)(radk) drivaa 
a fust order la«| sarvo. with a trnc cunatant of 
1 / 2 ^ aacon.1%. to <^nr:ata the actual rudr>rr poal- 
tion 5,(t)(tud>). Tna actual aliaror and ruddar 
puBltlor, f,(t) ar.d than asetta thr four 

‘natural* latatal dyr.arrca ftatr variables. na«a- 
ly roll-rate f(t)(rad sac), yaw-rata r(t)(rad/ 
sac), sidaclic angle r<t)(rad). ar.d oank an|la 
4(t*'rad). Ir. addition, a sind disturbanca stata 
van bla w(t). saa A^f^niis A. drivas tha aqua- 
tions in tha saaw way as tna iidaslip variabla. 


(I 




(t» . p*U) , B*(t) 


VMM 

4* ... I* 




:<t) ^ rc«t) 


4* 


4 * 

sa re 


dt 


nta fnllowinq smbIibm valuaa wara 'tsad 

Itaalw lataral acoalaration, ay„^,w0.2Sq*a 

4Va 

Naaiatai roll rata, p ■ » (a., -a*) 


Naaiaua sldasllp anqla, f 


«V* 


‘v/io^ 


■ji 


Thus, tha stata aquations for tha lataral 
dynaaics are character izad by a v-diawnslonal 
atata vector alt) with conponants 

(1.4) a'(t)»li(t) r(t) ?(t) :(t) «^(t) 6^(t) 


Naaianas ban)i anqla, ♦ *0 . 6 rad (>45*) 

Waalwun coa w andad aileron rata^l.C! rad/sac 
Maaianai coaasandad ruddar ratawl.23 rad/sac 


e (t) t It) w(t)i 

ac re 

and the overall lateral dynarics take the fora 

(3.5) a(t)»A »(t)eh u(t)-:- ^(t) 


where the zero rwar white r.cisc vector ^(t) qen- 
rrateti t).c wird disturbance ar.d coriansates for 
aodallinq errors, once rore the ratrices 
chanqe with fliqht conditions (2). (131 while 


(3.6) 


._[o 0 0 0 0 0 1 0 ol 
[o 000000 loj 


3.3 Tha Lateral Cost Functional 


The lateral performance indca used (after 
several iterations) vaiqntcd the following 
variables : 


See [13] for an eatenaive discuanion of how thia 
parforaance criterion was derivali S 3 j and a 3 j 
art obtained fro* the open loop ^ mat rices 

There ia no natural way of arriving at a 
alagrlifiad a»dal for the lateral dynamics, as was 
the case with the longitudinal dynamics. Hence 
the banli angle cannot be allainatad. Although a 
bank angle sensor was deeawd taidaslrable, the 
weak obaervability of the bank angle caused large 
atata aatlaiation errors, using iCalsMn filters, in 
tha bank angle and tha sideslip angle if a bank 
angle tensor was not Included. For these reasons, 
it was dacidsd to asiploy a bank angle sensor and 
to penalize bank angle deviation, because bank 
angles larger than 20* can introduce significant 
nonlinaaritias through trigonoaatric functions 
(11. 


o lateral acceleration, a (t) 
o roll rate, p(t) ^ 

o sideslip angle. ?(t) 
o bank angle, $ (t) 

VS 


(in g's) 

(in rad/sec) 
(in rad) 

(in rad) 


o cosssanded aileron rate, {^^(t) 
o coaenanded rudder rate, c^^(t) 


The lateral acceleration, a (t), is not a state 
variable. However, for siraYl perturbations from 
equilibriuiR flight, it can be expressed as a 
linear comfination of the lateral state variables 
and the trim angle of attack, as, by the following 
relation 

(3,7) a (t)-^ I(ki-a»)p(t) + (k)!+l)r(t)4kie(t) 
y 9 

♦k*^(t)-iki6^(t) J-« (t) 

idtere the constants ki,..., kj can be found from 
the lateral open loop ^ matrix, and change with 
the flight condition. 


The following structure of the quadratic 
perforaiance criterion was establishadi 


Once mure, the 10 can be solved. Notice that 
the use of the perforaiance criterion (3.B) results 
in a atata weighting matrix g, (non-diagonal) 
which changes with flight conation. 

3.4 Siamary of Results 

The above perforswnce criterion gave reason- 
able responses for a variety of initial conditions. 
Its main characteristic is to reduce any lateral 
accelerations (by forcing the aircraft to go in 
coordinated turns) and to null out bank angle 
errors in a slower manner. 

Once sure we observed a constant daaiping 
ratio ( .515 ) for all auparsonic conditions and 

a relatively constant daag>ing ratio ( ,625 ) for 

all B\i>aonic flight conditions. No additional 
tradeoff studies wars conducted by changing the 
waists in tha cost functional. 
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4. S*naof*, ralMn Ftltvrs and 
Disciftr LQC i a*>p«n»atora 

4.1 Introduction 

Tha dl 9 ltal ia 9 >Iaaiantatton of tha control 
ayatan raqutraa th« dlacrata-t aolution of tha 
LQC |<roblaa (i0|. Aj wa ahall aaa )n tha naat 
aactlui, tha IdtAC approach raquiraa tha conatruc- 
tion of a Lank ol LQC controlleia, each of which 
containa a rtiacrata Kalnan filter (whoaa raald- 
uala ara uaad in probal'illty calcuiationa and 
wtioar atata aatiaatra arc uaad to ijanarata the 
adaptiva control aiqnala). Hence, in thia aection 
we preaant an overview of the laauaa Involved In 
the daaiqn of the IJOC controllera baaed upon tha 
noiay aanaor ncaauranents . 

4.^ ITie Sampling Interval 

A saaplinq rate of 8 ixeaturemanta/aecond 
waa aatabliahed. Such a alow sai^linq rate waa 
aelected so aa to be able to carry out in real 
tine the multitude of real time operations re- 
quired by the KKAC method. 

4.1 Senaora and Noise Characteristics 

As explained in the introduction, the quide- 
lincs for desiqn excluced the use of air data 
sensors. Thus, ateaaurenents of altitude, speed, 
anqle of attack, and sideslip anqle were not 
available. After some preliminary investiqationa 
it was decided that sensors that depend on trim 
variables (elevator anqle and pitch attitude) 
should not be used so as to avoid eatlmatinq trim 
parameters. Table III lists the ••nsors and their 
accuracy characteristics that were used in thia 
study. We streaa that the sensors s«asure the 
true variables every 1/6 seconds in the presence 
of discrete-zerc mean white noise with the stan- 
dard deviations qiven in Table III. 

Finally, we remark that in this study «se 
assumed that all sensors were located at the C.G. 
of the '.ircraft. 

4.4 The Dcsiin of Kalman Filters 

P'>r each fllqht condition the steady-state 
discrete-time Kalman filter, with constant qains 
was calculated, for both the lonqitudinal 
and lateral dynamic models. The level of the 
plant tfhite noise associated with the wind dis- 
turbance qeneration was selected so that we 
assumed that the aircraft was flying in cumulus 
clouds. (See Appendix A.) 

The decision to use steady state constant 
gain Kalman filters was made so as to stinisiize 
the computer memory requirements. 

Finally, we remark that in view of the slow 
sampling rate, the continuous time filtering 
problem was carefully translated into the equiva- 
lent discrete prc^lem ( 13 ] to (IS]. 

The constant covariance Mtrices of the 
Kalsian filter residuals, denoted by Sj^ioN' SiLAT 


for the longitudinal models and lateral models 
more computed for each flight condition denoted 
by i. Aa we ahall aee thaee ara important In tha 
ganaratlon of tha WUC variablaa. 

4.5 Itta Oaeign of tha Discrete LCC Compenaators 

Through tha uae of the eeperation theorem one 
can deaign the discrete LOG cosgienaetors . Thia 
ii^liad that the iJp problem defined in rontinuoua 
tisM in Sections i and 1 had to be correctly 
transformed into the equivalent d:screte-tiM 
problem In view of 'he alow measurement rate, 
effectively, we have used the transformations 
given In referen < (111, (11) to (IS). 

4.6 Racepitulation 

For each flight condition, indexed by i, a 
completa dlacrcta-tisw, steady atata, LOG compen- 
sator was designed for both the longitudinal and 
lateral dynanica. Each c oeye n aator genarated 
avery 1/8 racond the optlisal control, namely the 
optisksl coainanded elevator rata c^-(t) for the 
longitudinal d^amics, and tha optlsul cci>rtandad 
aileron rate ^^^(tl and rudder rate o^^(t), based 
upon the noisy measurements of the appropriate 
sensors (Sae Tabla III) every 1/8 second. 

Because of the appropriate traj.sformations 
of the continuous tiam VQC problan to tha discreta 
one, we noted no significant degradation in per- 
fotmanca at this low sampling rata. 

The need for adaptive control is obvious l>e- 
causa if we assuM that tha aircraft is in flight 
condition i, but wa use tha LQC coopenaator ob- 
tained for flight condition j for feedback con- 
trol, this missMtehinq may qenarata either an 
unstable eyatem or, often, a system with degraded 
parfonsanca . 


S. The M1AC Hathod 

5.1 Introduction 

In this section we present the basic idea 
behind the MHAC method, and discuss how it was 
in the F-8C context. In particular, we 
deaonstrste how the information generated by the 
lateral and longitudinal sensors is blended to- 
gether. Finally we make some remarks associated 
with tbe MMAC swthod and its general applicability 
to the design of adaptive control systesm. 

5.2 Ti\e 3aeic Idea 

Sunpnae one has N linear, discrete-time 
stochastic time-invariant dynasiic systesm, indexed 
by i-1, 2, ..., N, generating discrete-time 
measurementf. corrupted by white noise 
Suppose that at t>0 "nature" selects one of these 
systestf and places it inside a "black box." The 
true syctem generates a discrete set of measure- 
ments s_(t) . The objective is to apply a control 
signal u (t) to the true axxSel. 


S 


wvraion of tha MNAr aaUtod M^loyod la 
M folluwai on* cor.> tiuct> a dl«cr«to*t IM staady 
•totr 1 ^ cofitiollai for aaLh itoalali thua, ofM- 
kaa a tank of M ^unfair datura . Aa ahoam in 
flfuia 1, aacfi LQi> cor|«naatot ta dr I van by tfa 
actual control a^i^liad to u.a syatan, uit) , and 
drivafi by tha actual noiay i^asuranant vactor, 
£ltt. Theta ara two atijnala of intarast that 
aach coa(>anaatci <|araiata» at tina t 

( 1 ) tha control vector ^It), which would ta 
tha of>tiral ccr.trol if indaad tha ayataai 
in tha blai k boa (vi<. aircraft) was 
idantiral to tha l~tl. nodal 

(2) the 'aaidoal or innovations vector 
r^lt) 9 enaratad by each laiaMn filter 
(which la inside the i*th 14 c coapenaator) 

It turns out that (see references (41, (S|, 
|b|, 17|, | 8 | for eaan|.le) that fron the residuals 
of the falsur. filters one can recursivaly generate 
M discieta tire sequences denoted fy K (t), lei, 

2 ..., N, t>0, 1 2, ..., vhich under iuitable 

assusFtions are the conditional ptef liti** *t 
tlSK t, qivcr. the past neasuranarts £(tl, t't 
and controls u(J> , d^t- 1 , that tha i-th nodal la 
tlia true one. 

Assuieinq than mat these probabilities are 
generated on-line (the forrula will be 9 iven 
later) and 9 iver. that each VQO c o a p ansator qan- 
arates the control vector ^(t), then as shown In 
riquie 1, the WtAC method coeputss tha adaptive 
control vector u(t), which drives t)ia true system 
(vie. aircraft) and each of tha Kalman filters 
inside t)i« I 4 C corper.sators , by probabilistically 
walqhtinq the controls (t) by the asaociatad 
prcbabilitias, l.a., 

N 

(5.1) u(t)-£ P (t)u^(t) 

1-1 

S.3 Calculation of the Frobabilitias P^(t) 

Me assumw that at t-0, i.a., bafora any 
swasurements Are obtained, one )ias a sat of prior 
probabilitias 

M 

(5.2) Pi (0), ..., P (0), P (0)>0. I F (0)-l 

‘ ■ 1-1 

that raprasent our ‘loast qvmss* of which modal ia 
indead t)ia trua ona. 

In our varsion of tha fWAC sw'iiod sia hava 
aval labia t*.a staady-stata (constant) covariance 
SMtrik Sj of t)ia residuals sssociatad with tha 
i-th Kalswn filter. These N residual covariance 
sutrices are preconputable. let r denote tha 
p.usber of sensors 1 than we can precoa^uta tha 
H scalars 

(5.3) Bj*-((2v)'dat 

Prom t)ia residual vactor r^(t) qanaratad by each 
KalBMh filter we qenerate on-line tha N scalars 

(5.4) SI, (t)Sr,'(t)S,-'r,(t) 

i —1 — 1 ^ 


Iham tha prebabi litlas at tis« t, , 1 - 1 , 2 , 

M ara computed recursivaly from tlia probabil- 
itlaa at time t-l, Pi(t*l). by tha forsaila 
P (»-l)i •aapf-m (t)/2) 

(5.5) r^U)-;p » 

I P (t-l)B /aiVif-m (t)/ 2 ) 

1-1 ’ > A 

with tha Iritlal probabilitias , P '0) qivan. It 
has bean ciaiamd that (5), ( 6 ), ( 8 ), under suit- 
able assim^tions that aaympf°f Ically tha true 
model la identified with probability 1 . 

5.4 Important Maawrka 

1) It has been s>»awn by Nillner ( 8 |, that 
the (MAC smthod, i.a., qanaratinq the control via 
(5.1) IS not optimsl (it is optlasl undsr suitAla 
assumptions for tha last staqa of tha dynamic 
proqrasmiinq alyorithm). 

2) Tha rathe alyoritlim is appealing in an ad- 
hoc way because of Its flaad struct urs snd bacaoss 
its rasl-tima and mamory raquiraswnts ara v-adlly 
coaqiutabla. 

3) In tha version us.'d in this study, becstsw 
ws use stssdy-stata Kalman filters, rsthar than 
tlsm-varyinq Kalman filters, ths 

P (t) ars not exactly the conditional 
ptobabllltias. 

4) Me hava bean unable to find in tha cited 
literature a rigorous proof of convergsnee of tha 
claim that i n daad tha probability aasociatad with 
tha true modal will asymptotically converge to 
unity. 

5) Prom a haurlatic point of view, tha ra- 
emrsiva probability foranila (5.5) ma)ias eanaa with 
raspact to Idamhificatlon. If t)ia system is sub- 
ject to soam sort of paralstant axitatlon, than 
ona would expect that tha raalduals of tha Kalman 
flltar aasociatad with tha correct modal, say t)ia 
i-th ona will be *ssiall,'* whlla tha residuals of 
tha mlsswtchad Kalsmn filters (VI, j-1, 2, .... 

■) will be ‘large." TTius, If 1 Indexes th# 
oorract model mm would axpact 

(5.6) m^(t)«m^(t) all VI 

If such a condition praalsts over several swasura- 
menta, tha analysis of (5.5) stiows tiiat tha 
‘corraci;" probability P|(t) will Increase while 
tha ‘mlMstchmd modal” probabilitias will dacrasxe. 
To sea vhis ona can rewrite tha formula (5.5) as 
follows 

(5.7) p, (t)-r.(t-i)-r I P,(t-l)B;a«V<-*4<‘>/2jl‘^ 

‘ ‘ L-1 1 1 ^ J • 

Pj(t-1> ^(l-P^(t-l))Pjaxp{-mj(t)/2lj 
- r P U-l)B;axp{-ai (t)/2) 

l»»i ’ ’ ^ 

Otodsr our sssisgitions 


6 


(».«) 


• a|>{ -•^(0/21X1 

IS.«) 

i 


H«nc« tiM coftact ptofcabi Itty will 9 row acconllnq 
to 

y (t-n ii-p u-iiP • 

<».io» — >o 

£ r (t-ll6 ••ap<-« (t)/2> 

,-i J ) 

which that at l'|(t)'*l. tho rot* of 

growth alow* lown. 

Or th* uth*( hand, tor th* incorroct aod*!*, 
ind*K*d ty )^i, th* taoic oaa^tiono yi*ld 

-P (t-l)P (t-l)8 • 

(5.11) P (t)-P (t-l)= -3 i < 0 

’ ’ N 

£ P (t-l)8j**p{-«^(t)/2} 
li-l 


»o thot 1 . 1 * f robohilitics d*cr*a**. 

The »oR« conclusirri .told if w* rewrite 
(5.7) It, th* f.»rp 

(5.12) P (t)-P (t-l).f I p (t-l)6**Kf{-oi (t)/2)|"^ 

jp (t-1) I P (t-l)(’ •*«?(-» (t)/2) 

IPX ’ W 

-«,**«Pl--,tt)/2|)J 


The above discussion points out that t)iia 'idanti- 
ftcation" sch*B* is crucially dependent upon th* 
reqularity of the residual behavior betwrer the 
'■tatched' and ‘rissiatched* Kalsian filters. 

(6) Th* 'identification* schesw. In tersis of 
the dynamic evolution of the residuals will not 
worli very well if for whatever reason (including 
errors in th* selectiot. of th* noise statistics) 
th* residuals of tho Kalstan filters dc not )iav* t)e 
above reqularity ar sumptions. Tobe specific, sup- 
pose that for a prolonged sequence of S Ms a ura sT ts 
th* Halman filter residuals turn C/ut to ba such 
that 

(5.13) Bi|(t).m 2 (t)w...«si^(t) 


Then 


(5.14) *xp(- «^(t)/2} ' ® for all 1 

Under these conditions and (5.12), we can see 

that 

P (t-1) £ P (t-1) (B *-6 •)a 

(5.15) P (t)-P (t-1). jPl ^ ^ 

* N 

£ P (t-l)B.*o 

J-1 ^ ’ 


P (t-l) £ (B •-•.•)P. (t-l) 

* 1^1 * ^ ^ 

H 

I p It-DB • 

!•! ’ * 

I sppese that it tarns oat that one of th* B.*'s, 
and to ba epecific B|,*, is doadnant, i.*., 

(5.16) ■»» lA 

In this casa, th* PMS of *q. (5.15) will be nega- 
tive for all ifk, which means that all the P (t) 
will decrease while th* probability P|| (associated 
with th* doninant B|,*) will increase. This be- 
havior is very important, especially when one lies 
to t)ie mathenatics , and it has not been diacuased 
previoualy in th* literature to the beet of our 
knowledge. 

5.5 Application to t)»* P-8C 

The MMAC method can be used In a straight 
forward manner using either th* longitudinal or 
lateral dynamics of th* P-0C aircraft since we 
liave designed both longitudinal and lateral ICC 
ooagienaators for the available flight conditions, 
aa w* ra m tr lied in Section 4. 


On t)i* otlter hand, w* obtain indapandant 
information from the longitudinal and lateral 
ayatasM for th* saam flight condition (i.e. , nod*!) 
indexed by i . Hence , it sltould be possible to 
blend thla coat-.ined information into a eat of sin- 
gle probabilities. 

Under th* assumption that ths longitudinal 
and Istsrsl dynamics era decouplsd )(-P Dunn de- 
rived th* follnwinq relation. 


£i ***^ £i lAT denote the resldusl 
covariance suitrices of th* Kalauin filters, for 
the i-th flight condition, associated with t)i* 
longitudinal and lateral dynasdes respectively. 
Define 


(5.17) B^*^- 


(2«) 


- 1/2 


(5.1.) B,*^^- 


^UiT 

aet*4 


- 1/2 


where end are tha. nueber of longitudinal 

and lateral sameors. let r^ A>xd U^T(t) 

denote th* icalmnn filter reaiduml vectora at tine 
t, for flight oomdition i, associated with th* 
longitudinal and lateral dynamica respectively. 
Define 


(5.19) m^ LOn“> 


(5.20) 


s ^r^ (t).— * r 

iLAT ^ LAT'^^'^i LAI*4 


’-AT 


(t) 


Then th* overall probability that th* aircraft is 
in flight condition i at tine t, is gansrated by 
th* recurs iv* formula 


7 


a. 




aik 






eoiitrlbwt* « o— und*r«t«n41fi4 uf>on th* NNAC awthcd 
M • dMi<in ronc*^ . 


lum' iLnr 


•**’* ‘ ^ '"lUkT 




it)/2} 


Tti* duaunar.c* di»cuk>«4 aiiov* now rtfars 

to the ralatiw M^nttoda of 


(5.22) e ‘‘S • r • 

' Pj UM 1 lAl 


Ofevioucly tn* Mthod should ha •■|'«ctad to tfork 
Mil wfiar. t>oth looqltudir.al and lataral ir^lrM 
filter* «rc correctly dasiqnad *o that th* rasid- 
ual* of tha ’■vatchad'* falnan filter* are knallar 
than those of the ‘‘mKrtatched* ones. 


5.6 Discussion 

It should be innadiataly obvious, that If the 
MKAC ncthod IS control cf the r~8C 

aircraft (cr any ether physical cysten for that 
Mtter) , onr violates a multitude of theoretical 
assiaaptions . Th* effect of these u(>on the perform* 
anca of the overall system is difficult to estab- 
lish on ar. analytical basis, because the MT*>C 
system, in spite of its simple structure, repre- 
sents an cxtreMly nonlinear system. Hence, on* 
hat to rely on extensive simulation results in 
order to be able to make a Judvnent of th* per- 
formances of the overall alqcrithn. 

Since the aircraft never coincides with the 
Biathcratical models (recall the discussion on the 
differences in the data giver, in references (21 
and (1], the P^(t) are not truly posterior proba- 
bilities. Father they should be inte .preted as 
tis<e sequences that have a reasonable physical 
interpretation. Hence, in our opinion, the eval- 
uation of the HHAC method solely by the detailed 
dynamic evolution of the P^(tl la wrong. Rather 
it should be judged by the overall perforsiance of 
tha control system. In the case of tha regulator, 
this is easy since ore can always coagiarc tha 
response of the HMT.J system with that which was 
designed explicitly for that flight condition and 
coagiare the results. 

Me resMrlr that such a coaiparlson is snich snra 
coeiplex when one attaclis the case of pilot inputs 
which result in several conssanded swnauvers. 

These aspects ar* still under investigation. 

There are several unresolved problems as yet 
which pertain to the total mjnber of ax>dels to be 
used at each instant of time, how these siodals are 
to be selected, how they should be scheduled in 
tha absence of any air data, and how one can arrive 
at a final desigr. th«it meets the speed-mesiory lim- 
itations of the IBM AP-101 computer wlilcik is usad 
in the VASA F-8C DFBW program. 

Me hope that some of the sianilation results 
and discussion presented in the sequel can 


u tlaulatiom Hesmlts 
k.l Imlroduction 

A variety of siaxilations Itave been don* using 
both a linear model and nonlinear sxidel of the 
r-K aircraft, lhase siaulatior.s results are ty- 
pical. They are selected such that they ran dao- 
Ofuitrata 

1) the speed of identification of th* 

MMAC algurithoi 

2) the overall perfurmanc* of th* NNAC 

systeoi and 

1) t)t* B* doaiinant behavior discussed in 
Section 5. 

Sow remarks about th* M1AC method ace given in 
th* conclusions. 

6.2 The Siaailatlon Results 

The simulations were conducted at a )>igh al- 
titude (4C.OOO ft), supersonic (Mach ! 4) flight 
condition (F/C 419 in Table I). No plant noise 
was introduced. All .■!x>dcls available In th* MMAC 
controller were given equal a priori probabilitias 
being th* true model. 

Emperismnt il i 

This It a set of llneer simulationa with two 
degree sideslip angle ( a B-gust) at tisw t>0. 

No sensor noise waa introduced and th* Kalmar, 
filtara wars set at t)t* correct Initial condition* 

rigur* 2 allows tha prol>ability changes while 
th* aat uf models availabla In tha HHAC controller 
were P/C B, 14, 18, 19 an^* 20. Not* that th* 
true flight costdltlon was included in tha con- 
troller. Th* correct modal is Initially chosen 
with high prohability within a very abort period 
of time (less than 1 smc.) and than switcltes to 
another axidal aiowiy after a few seconds. Lateral 
acceleration is reaxived within about one second, 
while roll rate and sideslip angle ar* reduced to 
aero alamt as fast. Mith no noise perturbing the 
■ysteai, the states of tha uystam have settled to 
near aero after about five second*. Thus th* 
residuals in all the mismatch stable filters 
approach sero. In this case t)ia B* doadnant be- 
havior dlacusaed in Section 5 occurred. Figure 3 
shows ihe probability changes when th* true model 
(F/C 119) was not Included (which was substltutsU 
)>y '.'/C tl7) . M* observed tha uam* B* dominant 
behavior after about five acconds. The sost ispor 
tent point to not* ia that rusponsas of tha ig(AC 
system ar* alsmnt idtuitical. Figure 4 sttowa t)ie 
responses of lateral acceleration with and withr>u'c 
F/C 119 in Uie controller, respectively. 

Similar results wars obtained with other ini- 
tial conditiosw. Howe w sx , tha spaed with which 
th* 6* starta to oaminet* varies greeily. For 
example, with a roll rat* initial condition, it 


*‘~i1 * HK-h »uOK»i . lhar* ta vary 

Itttla dv^adat i«Mt »n t.ta ovarall ayataa parfer- 


»aB<arta«:t.< »< i 

?:•># ta a aat of twtr.linaar aiiailaticna vith 
an initial aia da^'ea an<|la of attack lar i-^atl. 
Sansot noiaa vaa inttoducad aid tha kalnan flltart 
vara aat at taro initial condltiona. 


f;aura S ahowa tha ct>ar*ipr of pruLat.ilit lar 
vlian ttM- aat of aodala availaMa in tha IMAC con- 
trol lai varr r/C 14. 17, ly, anJ ^0 and vttar. tha 
trua flight condition (F 'C 14) via incladed. Tha 
prottahilitlaa ara nora activa tnan thoaa w hava 
aaan In Laimriawnt tl. It it taliavad that tha 
faat variation of thaa« yroC>aM ) it laa la daa to 
a roafcinatlon of tha tianaiar.t raaionaa of tha 
ayataa and tha noiaa ai'tjvirncaa on tha aanrora. 
Movavar, tha triM flight condition la idontlfia<^ 
in about 1 second. T'a argla of attack raturna 
to It a ttlBVnrd valua kiu.in about i aaccr.da, tdula 
pitch rate ar.l ixriwl acca Iciat lor. ara reduced to 
raro alatoat aa faat. Ir thia caaa tha f* doaiinant 
bahavior onl^ occura for a vary abort period of 
tine. Itecauaa of the aenaor noiaa, it la not 
certain if tha drifting of prc'babi litiaa ara 
aainly due to the dovirant i* . Figure f ahova 
tha probability changer; when F/C *19 (trua) waa 
aubatitutcd by F/C *16 ir. tha KKAC controllar. 
Again, f ic renponaea of the ryjyC. ayatan are 
alaioft identical. Figure 7 ahova tl« reaponsea 
of the angle of attack with and without F/C tl9 
In the controller, reapactively . 


iMMfltiadiMl dMi^ and Na. I. Svyall did Mch of 
klM pro^rrMlnf . 

•a or* loda>tod to Hr. Jorroll ■. Illlett of 
MIMIAJC tdte n. groot HBoltor previdod countlooo 
heora of toe**'^ diacuMiono and dlroctloo a«d 
hoXpod iM formiioco tho porforaaace crltorio. 
to addition, va oro Indafctad to tho folloviaf 
Btaff of HMA/IJC for thair holp. crltlelaH and 
aopporti R. nontgaaary, J. Cora. C. Maolay, 
k. Ichy and k. Evnna. 
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fcjiyn4l» A 


Hind Dit-.urbA'' « Hodvi 


A* r*Mik»d lA SActioftk 2 And i, a contimaouA 
tlaw wind ditturbAiM;* nodt l ma includAd tot both 
tlM lAtAfAl And lon<|itvMlir.Al dyrtAOicB. corr#A- 
»endin« to A AtAtA VAfiAblr w(t). In this Appon- 
AIB WA 9 IVA thA MthAAIAtlCAl dAtAllS Of thlB MOdAl 

wttlfh WAA kindly pfovidAd by Nr. J. Elliott of 
MAAA/LAl. aa a rAAAonAblr aj proAiantion to thA von 
RAxnAn nodAl And tliA tUirwA Appro aiMt ion. It la 
ti^iortAr.t to rAAliXA thAt thA wind diAturtiAncA 
■ndol chATt9AA frooi flight condition to flight 
cenditicm. Tho i«wwr ApoctrAl donaity of tho 
wind diAturbAi.cw it qivAr. hy 



whort L, thA acaIa Itnyth, it 

I 200 ft At aaa IavaI 

(A. 2 ) L • AltitudA > J 500 ft 

/ linoArly intcr(«lAtA<* ’ ' twoon 

Vg la thA tpAc-d of AircrAft in toe, w in 
rnd/aec, and 


1 6 ft/AAC nonBAl 

(A. 3) 9 ■ I 15 ft/AAC in cuBMlua clouda 

I 30 ft/AAC in thundAratoraw 


To obtAin A AtAtA variAblA BodAl, A nonaal- 
lAAd AtAtv vAriAblA w(t) (in rad) la uAAd aa tha 
wind AtAte for both lAtaral and longitiidlnAl 
dynaoica. Tha ata' • variAblc w(t) la tha output 
of A firat ordar /ataai dr Ivan by centinuous 
whlta nolAA ((ti with raro awan. Thus tha dynaa^* 
lea of tha wind diaturba:.cc aaidel ara 9 iva by 


(A. 4 ) w(t) 



C(t) 


whora C(t) la taro naan whita noiaa with unity 
covarianca function 

(A. 5) E |c(t)C(T)| - 6(t-T) 

Tha daaiqn waa obtainad for tha intaraodlata caaa 
O • 15 (cuMilua clouda) . 


for the longitudinal dynaadca tha wind atata 
w(t) Influancaa tha dynamica In tha aaaw nanner aa 
tha angle of attack. Thus, in tha longitudinal 
atata aquationa tha wind state w(t) antara tha 
aquationa aa followa 


q(t) - .... 

. A A . ♦ 

a w(t) 

• 


1 1 

w(t) - .... 

.... A 

a w(t) 

• 


II 

a(t) - .... 

.... A 

a w(t) 


whara a , a , a cah ba found froa tha open 
IS II II 

loap loMfitadinal A aatTfa (2). 

In tha lateral dynaoica tha wind atata w(t) 
inCl a aneaa tha dynaoica in tha aaao nannar aa tha 
aidaalip angla. Thua, in tha lateral atata 
aquatioAM tha wind atata w(t) enters tha aqua- 
tiona aa followa 

1 p(t) - A a w(t! 

. •• 

in.M r(t) - ♦ Aj,****) 

I i(t) « A a w(t) 

wtiara a , a , a can ba found trom tha open 

II 11 II 

loop lateral A natrin (21. 


